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The magnetic dipole moment (MDM) and the electric dipole moment (EDM) of leptons are 
calculated under the assumption of lepton flavor violation (LFV) induced by spin-1 unparticles with 
both vector and axial-vector couplings to leptons, including a CP-violating phase. The experimental 
limits on the muon MDM and LFV process such as the decay l~~ — > are then used to constrain 

the LFV couplings for particular values of the unparticle operator dimension du and the unparticle 
scale Au, assuming that LFV transitions between the tau and muon leptons are dominant. It 
Q . is found that the current experimental constraints favor a scenario with dominance of the vector 

■ couplings over the axial-vector couplings. We also obtain estimates for the EDMs of the electron 

' and the muon, which are well below the experimental values. 

<N 



> 



- 1—1 

x 



PACS numbers: 13.40.Em, 12.60.-i 



43 

I ■ I. INTRODUCTION 

Scale invariant quantum field theories cannot interpret matter in terms of particles. Motivated by the Banks and 
Zaks model Georgi [2|,|3( conjectured a scenario that, unlike the one posed by the SM and its extensions, introduces 
scale invariant stuff associated with fractionary anomalous dimension operators. Georgi suggested that a yet unseen 
scale invariant sector may exist in the high energy theory such that scale invariant stuff may interact weakly with 
the SM fields. Such a hidden sector may manifest itself at an energy scale A^ > 1 TeV, and since physical particles 
cannot exist in this sector, the interactions with the SM fields would occur through scale invariant fields known as 
unparticles. Although a comprehensive study of this class of theories is very complex, it is possible to describe its low 
energy effects through an effective field theory. This allows one to study the phenomenological effects of unparticle 
. stuff. 

The appropriate theoretical framework to describe unparticle physics is the one introduced by Banks and Zaks 
[l|. The hidden sector is a BZ sector and the associated fields are described by renormalizable Obz operators. It is 
assumed that these fields interact with the SM fields through the exchange of heavy particles at a very high energy 
M.U- Below this energy scale there is nonrenormalizable couplings between the fields of the BZ sector and the SM 
fields, 1 which generically can be written as OsmObz / Mif M+dl3z ~ 4 - Dimensional transmutation is caused by the 
renormalizable couplings of the BZ sector at an energy scale Ky as scale invariance emerges. Below this energy 
scale, an effective theory can be used to describe the interactions between the SM fields and the BZ fields, which are 
associated with unparticles. The effective Lagrangian can be written as [HQ : 

^dez — du 

ClA = C °" M dsM+d BZ -4°SMOu, (1) 

where Cq u stands for the coupling constant and the operator dimension du can be fractionary. From theoretical 
considerations Q, it has been noted [EEl that unitarity is guaranteed in the interval d u > 1. The Lorentz structure 
of unparticle operators is nontrivial but it can be constructed from the nature of the primary operator Obz and 
its transmutation. Such a Lorentz structure can be scalar, Ojj, vector, O u , spinor or tensor, O^ . The respective 
propagators of these unparticles along with their interactions with the SM particles have been studied with detail in 

UMm. 
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1 The dimension of the respective operators are dgz and dgM 
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Shortly afterwards the unparticle idea came out, the study of its phenomenology was eagerly addressed (71 llOL lllj. 
One interesting effect could arise from the interference between unparticle and SM amplitudes, such as could occur in 
the Drell-Yan process at the Tevatron and the LHC [1(3, EH 1 m particular, the most peculiar effects could be observed 
in the invariant dilepton invariant mass distribution near the Z-pole [H| • This class of interference effects could 
also be evident in diphoton production at the LHC [13| . As far as the direct production of unparticles is concerned, 
it has been studied through mono-photon, e~e + — > yU, and mono-Z production, e~e + — > ZU [ill ], whereas the 
production of an unparticle accompanied by a mono-jet was studied in [141 . Several decays of SM particles into 
unparticles have been examined: t — > blA Z — > f fU [l(J[ll|, H — > yU [l5j . and Z ■fU [16]. In addition, other 
topics on unparticle physics have been studied, such as the possible effects of unparticles on CP violation [l?], [U[ , 
lepton flavor violation (LFV) [l9l - [2l| . and lepton electromagnetic properties [22T424T ] . 

In order to impose constraints on unparticle parameters, several experimental data have been used. The e _ e + — > yhi 
process was studied to explain j9i/ production at LEP 11]. It was found that LEP data are consistent with the values 
Ay = 1.35 TeV for du — 2 and Ay = 660 TeV for du =1-4. Unparticle constraints have also been obtained from 
experimental data on cosmology and astrophysics [25l - [28l |. 

Apart from the tree-level effects, loop induced unparticle effects have studied in the literature [2ll - [23l . l29l]. The 
electron magnetic dipole moment via scalar and vector unparticles was obtained in 0, I22J. This study was later 
extended for the muon magnetic moment due to scalar unparticles with LFV couplings [231 ] . whereas the lepton 
electric dipole moment via scalar unparticles was studied in (24|. In addition, the calculation of the fcrmion dipole 
moments via fermion unparticles was presented in 0] . The study of loop induced decays mediated by unparticles has 
also been addressed, for instance the decays k — > ^7 [U [23| and Z — > Ulj (29|. In this work we are interested in 
calculating the spin-1 unparticle contribution to the magnetic dipole moment (MDM) and the electric dipole moment 
(EDM) of leptons in the most general case when there is LFV interactions. To our knowledge this calculation has not 
been presented in the literature. 

The rest of the work is organized as follows. In Sec. II we present an overview of unparticle operators. Section III 
is devoted to the results for the lepton electromagnetic vertex mediated by vector unparticles, while the numerical 
analysis is presented in Sec. IV. The conclusions and outlook are presented in Sec. V. 



II. UNPARTICLES INTERACTIONS WITH THE SM FIELDS 



The interactions between the SM particles and unparticles occur through the exchange of heavy fields of mass Aiu- 
Once those heavy fields are integrated out, the effective Lagrangian that describes the interactions between the SM 
particles and unparticles is obtained. This effective Lagrangian is composed by a tower of effective operators that 
can be constructed out of the SM fields by invoking the SUl{2) x Uy(1) gauge symmetry. For instance, the effective 
interactions of spin-0 and spin-1 unparticles with SM fermions are as follows ]l|, 0, 0, H> [I3] : 



C "s - ^MiOu + ^IrffiOu, (2) 

£u v = ^JaM + ^z-Ji^m- (3) 

where i and j stand for the family index and Aj = Co u A^f z /A4ff M+dBz ^ 4 stands for the associated coupling 
constants. These effective operators break scale invariance. 

Due to the invariant scale nature of unparticles, their propagators can be constructed by means of unitary cuts and 
the spectral decomposition formula. Therefore, the scalar unparticle propagator is given by: 

AHP 2 )= 9 A t" A -p 2 -ie) d «- 2 (4) 
where the A c i u function, which is introduced to normalize the spectral density [ll| . is given as follows: 



_ 16tt 2 V¥ T(du + |) 
du (2tt)^u T{d u - l)Y{2d u ) [0) 



As far as the vector unparticle is concerned, its propagator is given by 
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A£V) = Af(p 2 )^(p), (6) 

where 7r p "(p) is given explicitly as 

7r^(p = - 5 ^ + a^. (7) 

The form of this expression is due to the spin structure of this class of unparticles . Also p 2 ^ M 2 , which is a 
reflect of the unphysical nature of unparticles. When the unparticle field is taken as transverse, it turns out that 
PiiK^ip) — 0, which translates into the condition a = 1. However, in the context of a specific conformal invariance, 
a = 2(du — 2)/(du — 1) @- In the limit du — > 1 + the propagator Ap(p 2 ) turns out to be the propagator of a massless 
scalar particle, as expected: 

lim A F (p 2 ) = ^. (8) 

d u ^l+ P Z 



III. UNPARTICLE CONTRIBUTION TO ELECTRIC AND MAGNETIC DIPOLE MOMENTS OF 

LEPTONS 

Among the best measured particle observables, the muon MDM, a M , stand outs: it has been measured with an 
impressive accuracy of 0.54 ppm. The current world average, which is dominated by the measurements of the E281 
collaboration at BNL, is given by [30l ] 

a® xp - = 116592089 (63) x 10~ n (0.54ppm), (9) 

where the statistical and systematic errors, 0.46 ppm and 0.28 ppm, have been added in quadrature. 

The SM theoretical prediction is given by the sum of the QED, electroweak and the hadronic contributions. A 
recent review of these calculations is given in [3l|. While the QED and electroweak contributions to a M have been 
calculated with a great precision, the largest uncertainty arises from the hadronic contribution, which is still under 
revision. The theoretical SM prediction is 

a SM = H6591834 (48) x 10" 11 . (10) 

where the a(e~e + —¥ hadrons) data have been used to calculate the leading-order hadronic vacuum polarization 
contribution [32| . There is thus a discrepancy between the experimental and theoretical predictions larger than 3.6 
standard deviations [331 ] : 

Aa M = af- v - - af* = 255 (80) x 10" 11 . (11) 

As long as this disagreement is attributed entirely to new physics [34|, the allowed minimal an maximal limits for 
these contributions, with 95% C.L are Aa M = (255 =F 1.96 x 80) x 10 1 . It may be however that such a discrepancy 
will reduce to an acceptable level once the hadronic contribution is determined with best accuracy. 

Another well studied fermion electromagnetic property is the EDM, df, which can provide an excellent probe of 
new sources of CP-violation. In the SM, the Cabbibo-Kobayashi-Maskawa (CKM) mechanism cannot account for the 
amount of CP violation required to explain the baryogenesis asymmetry of the universe. It has been long known that 
the experimental observation of an EDM of fundamental particles will hint to new physics as the SM predictions are 
very small. For instance, the electron EDM is predicted to be negligibly small as it arises up to the three-loop level 
via the CKM phase. Other models such as supersymetry, multi Higgs and the left-right symmetric model predict 
much larger values for d e [35|. The current experimental limit on the electron EDM with 90% C.L. is (36| : 

| del < 1-6 x 10~ 27 ecm, (12) 

whereas the experimental limits on the positive and negative muon EDMs with 95% C.L. are [37}: 

|d+| < 2.1 x 10~ 19 ecm, (13) 

|d~| < 1.5 x 10~ 19 ecm. (14) 

We will calculate the contribution to the lepton MDM and EDM induced by a spin-1 unparticle with both vector 
and axial- vector LFV couplings. 
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A. Lepton dipole moments via spin-1 unparticle 




f(p') u f(p) 

FIG. 1: Feynman diagram for the lepton electromagnetic vertex induced by unparticles. 

The Feynman diagram contributing to the electromagnetic vertex is shown in Fig. [TJ From Eqs. ([2]) and Q 
we obtain the interactions of the spin-1 unparticle with both vector and axial-vector LFV couplings. We will use 
Feynman parameters for the calculation and the spectral form for the unparticle propagator will be considered: 

^ _A du f°° dm 2 (m 2 ) rf "- 2 
2tt Jq p 2 — m 2 + ie 



a f(p z ) = " , , (15) 



After the momentum space integration is worked out, we will proceed with the spectral integral. 

After some lengthy algebra we arrive at the following results. The MDM of the lepton j due to spin-1 unparticle 
can be written as 



E E W\ 2 Fj(mi,d u ), (16) 

J— V,A i—e,fi,T 



where the Fj functions can be written as 



167r 2 sm(ndu) \A 



2\ d u-i 



FA™M = ^ 2 "« [w) fji^du), (17) 



with s/Fi — mj/rrii and 

fv{z,d u ) = —V / dx (l-x) d "- 1 x 2 - d "(l-z 2 x) d "~ 3 ((3(3-x)-4d u ) + z(3(d u -l)x + d u -3) 
^ - du Jo v 

+ z 2 ((1 + xfd u + (x - 5)x - 2) - z 3 (3(d u - l)x + d u - 3) a;) , (18) 

also fA{z,du) = fv{—z,du)- We note that there is a flip in the sign of the vector and axial-vector couplings. 

For our analysis below, we will also need the contribution to a M from a spin-0 unparticle with both scalar and 
pseudoscalar LFV couplings. We obtain a similar result as that given by Eqs. (fT6| and (fT7|) . with J running over S 
and P, while the fs function is 

fs(z,d u ) = -z f 1 dx (l-xf u x 1 - d "(l~z 2 x) d "- 2 (l + zx), (19) 
J o 

and fp(z,du) = fs(—z,du). These results coincide with those presented in [23| and serve as a cross-check for our 
calculation method. 

We would also like to note that in the case of diagonal unparticle couplings, we obtain for the spin-0 and spin-1 
unparticle contributions to af. 
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.1, _ A du Y{2-d u )T{2d u -l) m)\ 4(2 - <W t 2 



a, 



l6Tr 2 sm(ird u )T(d u + 2) \ A2 



2(4, - 2)|A^ + -^-rfW - 3|A^| 2 + (2^ - 1)|A P | 2 I ,(20) 



where Aj = Aj £ stands for the diagonal unparticle couplings. This result agrees with the result presented in [22j for 
the electron MDM. 

As far as the lepton EDM is concerned, we will consider the contribution from both spin-0 and spin-1 unparticles. 
The EDM of fermion j, which can only arise if both Ay and A^ (Ac? and Ap) are nonzero and have an imaginary 
phase, is given by 

d " = E E Im(A^)G (J , K) K,4,), (21) 

{J,K) i=e,fi,r 

where {J,K) runs over (V,A) and (S,P). The Grj^) functions are defined as 

a / 2 \ du — 1 

G '^-""- i »'- ^l(^)„„ (i) w^.*) < 22 ' 

with 

9<y,A){z,du) = n ~j /'dx (l-x) d "- 1 2 ; 2 - d »(l-z 2 a ;) d »- 3 (4(l-z 2 a ;)(2-^) + (l-3x) 
^ — «w Jo v 

+ z 2 ((a; 2 -l)d w + (x-l)a; + 2)), (23) 



9 { s,p){z,d u ) = - dx{l-x) du x l - d «{l-z 2 x) d u-\ (24) 
Jo 

In the limit of a heavy internal lepton, TOj <C m„ the integration can be dealt with and we obtain 

<W"*.*<)- 647rsin2Mw)TOi (^J , (25) 

r- < a \ i d u ~ l)duA du fmf\ " -1 

G(s,p){mi,d u ) = — — , 2 — -T2- • (26) 



647rsin (Trdn)rrii \ A w 

These expressions can be useful for the muon and tau loop contributions to the electron EDM and the tau loop 
contribution to the muon EDM. 



IV. NUMERICAL ANALYSIS AND DISCUSSION 



Apart from the muon MDM, experimental limits on LFV process are known to be useful to constrain LFV couplings. 
LFV processes involving the muon are among the most constrained by the experiment. For instance, there are 
stringent constraints on LFV muon decays: BR(^ -> ej) < 2.4 x 10" 12 [H and BR(/j -> 3e) < 1.0 x 10~ 12 fill- 
Furthermore, the bound on the \i — > cy rate is expected to be improved by about one order of magnitude by the 
MEG experiment [40(. However, less stringent constraints exist for LFV r transitions: BR(r — > /17) < 4.4 x 10~ 8 
(4lj . BR(r — > 3e) < 3.6 x 10~ 8 [42], etc. Below we will analyze the constraints on LFV spin-1 unparticle couplings. 



A. Muon anomalous magnetic moment 



A comprehensive analysis of the spin-1 unparticle contribution to the muon MDM would require to deal with several 
free parameters: six coupling constants, the unparticle scale and the unparticle operator dimension. We will take 
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another approach instead and consider some particular scenarios of LFV together with the hypothesis that there is 
no large cancellation between different unparticle contributions. First of all, we will assume that there is a hierarchy 
in the LFV unparticle couplings as it is assumed in other models of LFV, i.e., |A ep | < |A eT | < lA^I <C A" for all the 
unparticle couplings. Since experimental data strongly constrains LFV between the muon and the electron, we will 
concentrate on the scenario where LFV transitions are dominated by the couplings of spin-1 unparticles with the tau 
and muon leptons. More specifically, we will focus on the three following scenarios: 

A- Vector and axial- vector LFV couplings of similar size. 

B- Vector-dominated LFV couplings. 

C- Axial-vector-dominated LFV couplings. 

To discuss these scenarios we first evaluate numerically each term, Fj{rrii, du), of the sum of Eq. (|16|) . The results 
are shown in Fig. [2] as functions of the dimension du and for Au = 1 TeV. We do not show the contributions 
due to the electron loop as we assumed that the LFV couplings involving the electron are subdominant. While 
the contributions from vector couplings are positive, the axial-vector contributions are negative. We also note that 
FA(rrii, du) — — Fv(rrii, du) and so the vector and axial-vector contributions can largely cancel out. This effect is 
more evident for the contributions of the tau lcpton. Furthermore, the contributions from axial-vector couplings are 
larger in magnitude than the vector contributions, with the largest contribution arising from the muon loop. It means 
that as long the vector and axial-vector unparticle couplings were of similar order of magnitude (scenario A) the 
total contribution to the muon MDM could be negative. This is a scenario disfavored by the current experimental 
data, which requires a positive contribution to from new physics. The situation is rather similar for the scalar 
and pseudoscalar contributions, which for comparison purposes are shown in Fig. [3J although in this case the scalar 
contributions are slightly larger in magnitude than the axial- vector contributions. If all the unparticle couplings were 
of similar size, the unparticle contribution to a^f could be negative as the axial- vector plus pseudoscalar contributions 
could be dominant. However, even if all the unparticle couplings were large, the total contribution to a M could still 
be positive by a sort of fine-tuning. We consider the scenario in which the axial-vector couplings are negligible in 
comparison to the vector couplings (scenario B), such that the total contribution to is positive. Below we will 
analyze the possible constraints on vector unparticle couplings under this assumption. 
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FIG. 2: Partial contributions from the indicated lepton loop to the muon MDM due to vector (V) or axial- vector (A) unparticle 
couplings as a function of the du dimension and for Au = 1 TeV. The absolute values of the (negative) (A) contributions are 
shown. The horizontal lines are the minimal and maximal allowed limits on new physics contributions to a M with 95% C.L. 
Notice that the contribution from the tau loop are almost indistinguishable. 

If the current discrepancy between the theoretical SM prediction and the experimental value of the muon MDM is 
assumed to be due entirely to unparticle interactions, the allowed limits with 95 % C.L. for this class of contributions 
are 98.2 x 1CT 11 < afj < 411. 8 x 10~ n . We have found the allowed area on the \X^\ vs \X^\ plane, which is shown 
in Fig. [H for several values of Au and du that are consistent with the bounds obtained for the scale Au from mono- 
photon production at LEP [ll| (see Table H]). As inferred from Fig. [4j the experimental data allow a magnitude of the 
vector couplings as large as unity for du — 2 and Au — I TeV, whereas the most strong constraints are obtained for 
du close to unity. In general the limits on | X V T | are slightly stronger than the limits on | Xy^ | , which is in agreement 
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FIG. 3: The same as in Fig. [2] but for the scalar and pseudoscalar contributions to a M . The absolute values of the (negative) 
(P) contributions are shown. 



with our assumption. We also note that the bounds on the vector and axial- vector unparticle couplings are of similar 
order of magnitude than the bounds on the scalar and pseudoscalar unparticle couplings, as shown in Fig. [5] To 
obtain that plot we have assumed that the dominant contribution to the muon MDM arise from a spin-0 unparticle. 



d u 


Au [TeV] 


2.0 


1.35 


1.8 


4 


1.6 


23 


1.4 


660 



TABLE I: Constraints on Au from mono-photon production data at LEP with 95%C.L. and assuming Ay = 1 [ill ]. 



B. Bounds from the decay r — ¥ 3/Lt 



It has been long known that LFV couplings can also be constrained from the experimental bounds on the tree-level 
induced decays Ij — > ^jtj^t^. We will examine the bounds obtained from the r — > 3/i decay, which involves the Ay 1 
and Ay 1 couplings. The calculation for the l~[ — > tjt~^t~^ decay width was already presented in [19( for the scalar 
unparticle contribution. We have calculated the contribution from vector unparticles and the result is presented in 
the Appendix |XJ Let us consider the scenario we are working in, and neglect the axial- vector contributions. With 
these assumptions we can write for the r — > 3/i branching ratio: 



BR(r -> 3/i) 



2 8 7 



A 



du 



sin(d7r) 



m 7 



4(du-l) 



|AHW|V 



< 10" 



(27) 



with T r the tau mean life. The right-hand side of the inequality is the experimental constraint [33| and the r/i function 
is defined in Appendix [A"l 

Eqs. (fTfJ)) and (|27p can serve to further constrain the allowed values of the coupling constants. We show in Fig. [5] 
the allowed area obtained after numerical evaluation of Eq. (|27p for several values of An and du- We observe that 
the t — y 3/i decay constrains considerably the size of the Ay" parameter, whose allowed value is extremely small for 
du close to unity, where the allowed values are considerably larger. There is also an allowed area in which the 
opposite is true (|Ay"| 3> | Ay/ 1 1 ) , but we will not consider it as conflicts with our previous assumptions. For comparison 



purposes, we also have obtained the allowed area on the |A^ | vs |A^ | plane when it is considered that LFV scalar 
unparticle couplings give the main contributions to the muon MDM and the r — > 3/i decay. The results are shown 
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FIG. 4: Allowed area on the |Ay T | vs |Ay*| plane consistent with the experimental limit on the muon MDM with 95% C.L. We 
used the indicated values of Au and du- 



du=1.2, A=100 TeV du=1.4, A=10 TeV 
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FIG. 5: The same as in Fig. [4] but for the bounds on the |Ag M | vs. |Ap'| plane. 



in Fig. [7] for several values of and dy. We observe that the allowed size of the scalar unparticle couplings is of 
similar order of magnitude than for the vector unparticle couplings. 

We also analyzed the possible constraints obtained from the r — > fij decay. It involves the Ay T , Ay M and Ay 
couplings. We find that the respective constraints are less stringent than the ones obtained in Fig. [6l so we will not 
consider this decay mode. 
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FIG. 6: Allowed area on the |A^ T | vs |A^ M | plane consistent with the experimental constraints on the muon MDM and the 
t — y 3fi decay for the values of Au and du indicated in each plot, when the dominant contribution to these observables arise 
from the vector unparticle couplings.. 
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FIG. 7: The same as in Fig. |6j but when the dominant contribution to a M and r — > 3/i arise from the scalar unparticle couplings. 



C. Electron and muon electric dipole moment 



In order to have a nonzero EDM, both kind of unparticle couplings, vector-axial or scalar-pseudoscalar, must be 
nonzero. Also, an imaginary phase is required. We will try to get an estimate of the order of magnitude for the 
electron and muon EDMs induced by unparticles. Apart from the magnitude of the (Ay^j and |A^ p | parameters, 
the EDM depends on the associated CP-violating phases, so its analysis is even more complicated. The most general 
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form for the unparticle couplings is A 1 / = \X l j | exp(— i9 % j ), where % j = —O^J is a CP-violating phase. Therefore 
Im(AjA^) = | A j ||A^-| sin 58^ JK y with SO'^j K ^ = 8 l j — 8^ the relative phase between J and A" couplings. Depending 
on the relative sign of the CP violating phases, 89^ v A ^ and S9^ s p y the partial contributions can add coherently or 
destructively. We will content ourselves with considering a somewhat restrictive scenario, which will allows us to get 
an estimate of the order of magnitude of the EDM. First of all, as was the case for the muon MDM, we will assume 
that there is no large cancellation between different contributions to the EDM, so we can analyze each contribution 
separately. Also, to be consistent with the previous discussion, we will assume the following hierarchy |Ay T s | ^ |A^ e s | 
and \\yg\ ^ lAy 1 ^. It means that the EDM would be entirely driven by the tau loop contributions. In view of this 
scenario, we numerically evaluate the tau loop contributions to the G(j_k) functions of Eq. (f2l"j) . The results are shown 
in Fig. [5] We observe that the contribution from the tau loop to the electron and muon EDMs are indistinguishable, 
which means that the EDM is entirely controlled by the magnitude of the couplings constants and the imaginary 
phases. For this calculation we used the approximate formulas given by Eqs. ([25]) and (|26|) . and made a cross-check 
with the exact calculation obtained by numerical integration of Eqs. (|23|) and (|24|) . 




1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

du 

FIG. 8: Partial contribution from the tau lepton loop to the electron and muon EDM due to vector-axial (V, A) and scalar- 
pseudoscalar (5*, P) unparticle couplings as a function of the du dimension and for An — 1 TeV. The absolute value of the 
(negative) (V, A) contributions is shown. The curves for the tau contribution to the electron and muon EDMs are indistin- 
guishable. 

We now consider the constraints on the coupling constants obtained above to estimate the muon EDM in our 
working scenario. Both (V,A) and (S,P) contributions are of the order of 10 -17 e-cm for du ~ 1.1 and 10~ 21 e-cm 
for du — 1-6. We will consider the case in which the (V, A) coupling prevails over the (S,P) couplings. The allowed 
magnitude of the |Ay T | parameter is of the order of 10 -5 for du = 1-1 and 10 -3 for du = 1.6. A good assumption for 
the |A^ T | coupling is that its magnitude is about one or two orders of magnitude below \X V T \. This would not spoil 
the bounds we obtained above. Therefore the muon EDM is about \d^\ ~ | sin 5^ | x 10 -30 e-cm for du — 1-1, 

whereas \d^\ ~ | sin5^^| x 10~ 29 e-cm for du — 1.6. These results are well below from the experimental limit on 
the muon EDM as the imaginary phase is expected to be very small. The magnitude of the coupling constants for 
LFV transitions involving the electron are expected to be more suppressed as long as the hierarchy discussed above is 
respected. Then we can estimate that the electron EDM induced by vector unparticles is well below the 10 -30 e-cm 
level and far from the experimental measurement. Moreover, since the constraints on LFV spin-0 and spin-1 unparticle 
couplings are similar, the contributions to the muon and electron EDMs from both spin-0 and spin-1 unparticles are 
expected to be of similar size. Even if all different unparticle contributions add coherently, it is hard to expect a large 
size of the electron and muon EDMs due to unparticles. 

V. CONCLUDING REMARKS 

The lepton magnetic and electric dipole moments were calculated in the framework of unparticle physics assuming 
interactions mainly driven by a spin-1 unparticle with both vector and axial-vector couplings to leptons. Analytical 
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expressions were found that agree with previous calculations for the muon MDM in the limit of flavor conserving 
couplings. 

The muon MDM was numerically evaluated to obtain bounds on the coupling constants from the experimental 
measurements assuming that LFV is mainly dominated by LFV couplings involving the muon and tau leptons. It 
was found that the latest experimental data for the muon MDM favor a scenario in which the spin-1 unparticle 
contribution is dominated by the vector couplings since the contribution from axial-vector couplings is negative. It is 
also found that the decay r — > 3fi can constrain considerably the allowed magnitude of the Ay M and Xy couplings. 

As far as the EDM is concerned, since it depends on several free parameters, we content ourselves with estimating 
its order of magnitude. It is found that that the unparticle contributions are well below the experimental limits of 
the electron and muon EDMs. 
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Appendix A: Decay £r ->• l^t~^t\ 

In this appendix we present the calculation for the i~ — > £j£ k £ k decay mediated by spin-1 unparticles, which 
proceeds through two tree-level Feynman diagrams (in the first diagram the unparticle couples to the final £ k £ k 
leptons, while the second diagram is obtained through the exchange of £j £^). After introducing the Feynman 
rules for the spin-1 unparticle, the decay width can be written as 



2 /™ \ 4 ( d «- 1 ) 

m -+ ejr k £+) = 



2 9 7T 3 



sin(<i7r) 



\ a£ j J dxi I dX2 C^ 1 ' 2 + |Al2 ' 2 " 2Rc(All2) ) • (A1) 



If .7 = k a factor of 1/(2!) must be included as there are two identical particles in the final state. The integration area 
on the xi vs. x 2 plane is 

4sfe - s 2 < x\ < 1 - 2sj, 



> 1 / ((l^) 2 -^ 2 ) (x 1 + s 2 -As 2 ) 

X2 < 2\ 1 - xiT y xTTs* ) (A2) 



where Sj = rrij/mi. In addition, Mi arises from the squared amplitude corresponding to each Feynman diagram and 
M12 from their interference. These terms depend on the du dimension and the x\, x 2 variables, and are given by: 

l^il 2 = |A^| 2 (|A^| 2 /i(^,^) + lA^I 2 /^-^-,^)) + |A5f | 2 (|A^| 2 / 2 ( Sj , Sfe ) + |A^"| 2 / 2 (- Sj , Sfe )) , (A3) 

\M 2 \ 2 = |AVT (\\$\ 2 gi( S3 ,s k ) + \\^ gi (-s 3 , Sk )) + \X%\ 2 {\X$\*g 2 { Sj ,s k ) + |A^ fe | 2 52 (- Sj , Sfe )) 

+ 2Re((A^AV fc A^A^)).9 3 ( Sj , Sfc ) (A4) 

M 12 = \% (xf (x^X^h^Sk) + XfX^fni-sj, 8 k j) + X k v k (x^ k \fh 2 { Sjl s k ) + X*fX^h 2 (- Sj , Sfe ))) 

+ (V<r>A,h 1 (s j ,Sk)<r>h 2 {-s j ,-s k )), (A5) 
For the sake of clarity we omitted the explicit dependence on x\, x 2 and du- The gi and hi functions are 

h{u, v) = - (2u 3 + (l + Xl )u 2 + 2 (2v 2 + Xl )u+ ( Xl - 1) (2v 2 + x x ) + 2 (a* + x 2 - 1) x 2 ) (u 2 + Xl ) 2{du ~ 2) , (A6) 
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h(u,v) = f 1 (u,-v) + 2v 2 (2u + x 1 -l)(u(3u + 2) + 2x 1 + l)(u 2 + x x ) 



2d u -5 



(A7) 



9i(u,v) = - 



{v - 1) (x 2 (v + 3) + 2v 2 + v + l) + u 2 ( (5 - x 2 ) v 4 - ((x 2 - 10) x 2 + 4xi (x 2 + 1) - 3) v 2 

2x\v - 2v 5 + x 2 (-2x 2 2 + x 2 - 4xi (x 2 + 1) + l) ) - 2x\ (2v (v + 1) + 2x x - 1) - 4 (x x - 1) 2 

2uw (2w + x 2 - 1) (w 2 + x 2 ) [y [3v + 2) + 2x 2 + 1) v 4 - x\ ((4x x - 3) v 2 + 3w 4 + lOv 3 + 4 (xi - 1) Xi) 



- x 2 v 2 (6v 3 + v 2 + 4xi (2xi - 3) + 3) - 2x^ (v 2 + x 2 ) 



v2(du-3) 



(A8) 



92(«,«) = 9i(~u, ~v), 



(A9) 



ff 3 («, v) = -4 (u 2 - -y 2 ) (1 + x 2 ) (v 2 + x 2 ) 



2(d u -2) 



(A10) 



h\(u,v) 



i 6 (v + l) + u 5 (l - x 2 (v - 1) + v 2 ) +u 4 ((xi + x 2 - 4)v 2 + (xi - x 2 )v 

- Av 3 + 2 (x 2 + xi (x 2 + 1)) ) - u 3 ( - x 2 (v + 1) - x 2 (w (2(v-2)v-x 1 + l) + 3xi) 

+ w (w (v (4w + xi - 1) - 2xi + 4) + 1) + x 2 - xi) + u 2 (-2(x 1 + l)v 4 + (-3xi - 2x 2 + 3) v 3 

+ (3x 2 + (x 2 - 7) x\ + (x 2 - 7) x 2 + l) v 2 - x 2 (xi + x 2 - 1) v + 2 (x 2 - 1) x 2 + 2x x x 2 (3x 2 - 1) 

+ x\ (4x 2 + 1) ) - u( - xix 2 (v + 1) - xix 2 (v (2 (u - 2) u + 1) + 2xi - 1) 

+ u (v (xi (w (4u - 3) + 2) - u - 2x^ - l) + xi) ) - (x 2 - l) v 4 + (xi (1 - 2x 2 ) - 1) v 3 

- xi (x 2 - 1) x 2 v + (2x 2 + xi (2x\ + (2x 2 - 3) Xi + (x 2 - 3) x 2 + l)) v 2 



2x\x 2 (xi + x 2 - 1) (xi + x 2 ) 



K +Xi 



x 2 )' 



(All) 



/i 2 (u, i>) = - [w 3 (x 2 + 1) (u - 1) + m 2 ((xi + 2x 2 - 2) v 2 - (v 2 - l)v + x l + 2x 2 (xi + x 2 )) 

+ u ((2xi + 3x 2 - 1) v 3 - 2 (x x + 3x 2 - 1) v 2 + (x 2 (2xi + 3x 2 - 2) - 1) v - 6v 4 - x 2 (2xj + x 2 - 1)) 
+ u 4 {v + l)+xl(3(v-l)v + 4xi - 2) + (xi - 1) v 2 (v (3v - 2) + 2xi - 1) 



+ x 2 ((6xi - 4) v 2 + (1 - 2 Xl )v -3v 3 + 2 (xi - 1) Xi) + 2x 3 (u 2 + x x ) du 2 (v 2 + x 2 ) 



\du-3 



(A12) 



1. Decay li -> 3lj 

In the case of this decay, the above expressions simplify considerably. We can write the decay width as 

r(£, -> 3ij) 



with 



mi 






2 W TT 3 


sin(c?7r) 


UJ L 



+ 



|A#'l 2 (|A«|V( Sj ,^) + |A«|V(-*i,<fc)) 
|Ajl 2 (|A«| 2 tj 2 (* i) ^) + |A^ 2 % (- Sj -,d^))+2Re(Aj'A^ A« A*') (A13) 



r)i(z,d u )= J dxi J dx 2 £i(z,d u ) 



(A14) 
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where the integration limits on the x\ vs. X2 plane are now given by 

3s? < xi < 1 - 2s j, 



> 1 / /((I " xi) 2 ~ 4s?) (a;i-3«3)\ 

»l 5 ('-« \/- — 4rrf -J ■ (A15) 

and the £j functions, whose explicit dependence on x\ and 22 has also been omitted for clarity, are given by 

£i(z,du) = -\(2q 2 i du ~ 2) ((z + l){3z-l)x 1 +x 2 1 +2(x 1 +x 2 -l)x 2 + (6z-l)z 2 ) 

+ (<7i<?2) d "~ 2 (xi + x 2 + z - 1) (x a +x 2 + (3z + l)z) ) + (xi <-» x 2 ), (A16) 

6(^, d«) = 2 ( 4 9i d " _5 + I) 2 (si + 2z - 1) z 2 + 2^ rf «- 4 (( z + \f Xl -2x\-x\-2 {x x - 1) x 2 - (lOz + 3)z 2 ) 
+ 2qf u - 3 q$ u - 2 (z + 1) {xx + 2z - 1) (x x + 2z 2 + z) z 
- (<M2) d "~ 3 2 2 (z + l) 2 (xix 2 + Ol + ^)^ 2 + z (z (z + 1) - 1)) 

- (gi92) d "" 2 ((xi + x 2 + 2) z 2 + x\ + x\ + 2x lX2 -X1-X2 + (7z 2 - l)z) )+{xn> x 2 ), (A17) 

6(z,d w ) = -2(^"- 3 g 2 d "- 2 (l + x 2 -2z 2 )z 2 +^"- 2 g 2 d "- 2 (x 1 +x 2 -l)(x 1 +x 2 )) +(i x oi 2 ). (A18) 
We also introduced the notation q, ; = SEj + z 2 . 
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